Further evidence for linearly-dispersive Cooper pairs 
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Abstract 



A recent Bose-Einstein condensation (BEC) model of several cuprate superconductors is based on bosonic 
Cooper pairs (CPs) moving in 3D with a quadratic energy-momentum (dispersion) relation. The 3D BEC 
condensate-fraction vs. temperature formula poorly fits penetration-depth data for two cuprates in the 
range 1/2 < T/T c < 1 where T c is the BEC transition temperature. We show how these fits are dramatically 
improved assuming cuprates to be quasi-2D, and how equally good fits obtain for conventional 3D and 
quasi-lD nanotube superconducting data, provided the correct linear CP dispersion is assumed in BEC at 
their assumed corresponding dimensionalities. This is offered as additional concrete empirical evidence for 
linearly-dispersive pairs in another recent BEC scenario of superconductors within which a BCS condensate 
turns out to be a very special case. 

PACS numbers: 73.22.-f: 74.20.Mn; 74.72.-h; 74.78.-w 



Introduction 



A Bose-Einstein condensation (BEC) model was applied by Rosencwaig in Ref. pQ to address seven cuprate 
superconductors (SCs) with transition temperatures T c at optimal doping ranging from 22 K to 133 K. These 
are: La 2 _ K Sr :E Cu04 (LSCO), Nd 2 - K Ce :E Cu04 (NCCO), YBa 2 Cu 3 07-<, Y123, Bi 2 Sr 2 CaCu 2 8 _ a Bi2212, 
Bi 2 Sr 2 Ca 2 Cu 3 Oio- y (Bi2223), HgBa 2 CaCu 2 7 - y (Hgl212) and HgBa 2 Ca 2 Cu 3 9 _ y (Hgl223). His starting 
point is the well-known fact that BEC in an ideal Bose gas occurs below temperatures T such that the 
thermal wavelength A = hf ' \J2-Km,B^BT becomes larger than the average interbosonic separation, with tub 
the boson mass, and h, ks the Planck and Boltzmann constants, respectively. More exactly, BEC sets in 
whenever 

n B X 3 > 2.612 (1) 

where ns is the boson number density, and A is taken as the bosonic quasiparticle diameter. This leads to 
a critical temperature T c given by the familiar formula 

2irh 2 n 2 ^ 3 _ 3.31fi 2 n| /3 2/3 , . 

c ~ (2.612)2/3 ms fc B " mB k B B ( ) 

of conventional BEC theory. He identifies an interaction distance withA, which thus becomes T- independent, 
while the number-density tib of weakly-interacting "preformed" electron or hole pairs acquires a T-dependence. 
Associated with ([2]) is the BE condensate fraction 

where Nq(T) is the number of bosonic pairs at temperature T in the lowest-energy state with (total, or) 
center-of-mass momentum wavenumber K = 0, and A/q(0) is that same number at T = 0. 
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Cooper pair dispersion 



All of this assumes three dimensions (3D) and that boson excitation energies are given by 

s K = h 2 K 2 /2m B (4) 

with as would hold if the composite bosons moved in vacuo such as, say, a deuteron of mass tub — m p + m n 
in empty space with m p and m n being the proton and neutron masses. However, in the presence of the Fermi 
sea of the other single charge carriers the bosonic "dispersion relation" becomes linear [5] [3] in leading order 
rather than quadratic as in (HJ). Ref. [2] first mentioned, and Refs. |3][Zj later discussed this linearity in 
detail. It is associated with the original Cooper pair (CP) problem [8] of two electrons (or holes) above (or 
below) the Fermi surface of the remaining system electrons. It was also found in a more general view of CPs 
in Refs. 9 10J within a many-body Green's function formalism treating both electron- and hole-pairs on an 
equal footing. For either ordinary or generalized CPs, the leading term in the if -expansion is linear. This 
linearly-dispersive "moving CP" object is often confused in the literature with the more common Anderson- 
Bogoliubov-Higgs (ABH) [11] (Ref. [H] P- 44) [13] collective excitation which is also linear in leading order, 
but which is just the sound mode of the many-fermion system. By contrast, in a many-boson system these 
two modes, the "particle" and "sound" modes, are apparently identical [14] [15j. 



A particularly clear example comparing linear and quadratic dispersion is perhaps the analytical result of 
Ref. [B] in 2D for an attractive delta potential assumed between electrons. This interfermion interaction 
mimics the net effect of Coulomb repulsion plus attractive, say, electron-phonon interactions. The 2D delta 
potential well, which otherwise supports an infinite number of bound levels, is imagined "regularized" [T|}] to 
support a single bound level of energy — B^ as occurs [17], e.g., with the two-parameter Cooper/BCS [8] [19] 
model interelectronic interaction. Miyake used this interaction to obtain [18] both the zero-temperature BCS 
gap A and the chemical potential fi analytically in terms of £2. Since the regularized delta well turns out 
to be infinitesimally weak, its + strength can be eliminated [6] in favor of B2 which then plays the role 
of coupling constant, with < B2 < 00 spanning weak to strong coupling. Instead of Q, a more general 
analytical expression found in Ref. [5] that includes the Fermi sea is 



2 

ek = — fivpK 

7T 




h 2 K 2 

2(2m e ) 



0(K 3 ) (5) 



where vf is the Fermi velocity defined through the Fermi surface energy Ep = rn e v F /2 and m e is the effective 
electron mass. The leading term is linear, and only in the vacuum limit (vp — > 0, implying Ep — > 0) does it 
precisely become the quadratic ([4]) with mg = 2m e expected physically for any fixed coupling B2. Fig. 1 of 
Ref. [6] exhibits the smooth crossover in 2D from a purely linear to a purely quadratic form, as one increases 
coupling and/or as one "switches off" the Fermi sea medium (nonzero Ep) in which the pair propagates 
down to the pure vacuum (zero Ep) medium. A very similar behavior was also observed in 3D [7], but only 
numerically. 



Bose-Einstein condensation 



Expressions more general than ([2]) and ([3]) but reducing to them, are known [20] in any dimensionality d > 
(integer or not) and for any dispersion relation 



SK = C s K s 



with s > 0. 



They are 



T r = 



sT{d/2){2ir) d n B 



-, s/d 



2^T{d/s)g d/s {l) 



s/d 



(6) 
(7) 
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and 



N (T)/N o (0) = 



1 - (T/T c ) d / S for T <T C 



(8) 



otherwise. 



Here are the Bose integrals which for a > 1 coincide with the Riemann Zeta-function C(<r) and diverge 

for a < 1, r(cr) is the gamma function, and = N/L d is the d-dimensional boson number density. The 
divergence of <?o-(l) for <j < 1 ensures from that T c = for all d < s, but that otherwise T c is nonvanishing. 
In d = 3 and quadratic dispersion s = 2 and, if C2 = h 2 /2rriB, © and © respectively become © and ([3]), as 
53/2(1) = C(3/2) ~ 2.612. However, for s = 1 BEC can occur for all d > 1. This coincides, fortuitously, with 
all dimensions where actual superconductors have been found to exist, down to the quasi-one-dimensional 
organics [21] [23] consisting of parallel chains of molecules. As regards dimensionality, therefore, the BEC 
picture contrasts sharply with the BCS scheme where T c is nonvanishing for all d > even though no 
exactly ID superconductors have been found to date. In fact, beautiful experiments |24] )25| with nanowires 
of different thicknesses sputter-coated with an amorphous superconductor (T c ~ b.5K) have shown how 
superconductivity is extinguished for the smallest nanowire diameters interpreted as approaching precisely 



Although the creation/annihilation operators of BCS pairs do not obey the usual Bose commutation rules 
[see Eqs. (2.11) to (2.13) of Ref. [Hj; see also p. 38 of Ref. [2]], CPs in fact satisfy BE statistics. Indeed, 
BCS pairs and CPs are distinct. A BCS pair is defined with fixed total (or center-of-mass) momentum 
wavevector K = ki +k2 and fixed relative-momentum wavevector k = (ki — k2)/2, whereas a CP is defined 
with fixed K only, since a sum over k is implied in any conceivable formulation of CPs. This is because 
in the thermodynamic limit an indefinitely large number of BCS pairs, each with fixed momenta Tiki and 
h]i2, correspond to different relative momenta fi.k but whose ki and k2 add vectorially to the same total K. 
These remarks apply even when only K = pairs were considered in Ref. [19] 



Empirical evidence for the linearly- dispersive nature of CPs in BSCCO has been argued by Wilson [26] 

as being suggested by the scanning tunneling microscope conductance scattering data obtained by Davis 
and coworkers [27] [28) in this cuprate. In Figure 1 we present additional evidence, based on experimental 
data from penetration-depth measurements in two 3D SCs |29) [30] and two quasi-2D cuprates [31] [33], as 
well as from gap measurements in a quasi-lD nanotube SC [33]. When plotted as a presumably universal 
"normalized order parameter" the data depart substantially (at least in 3D and 2D) from the BCS normalized 
gap order parameter, but are seen to agree quite well, at least for T > 0.5T C , with the pure-phase (only 
2e- or 2h-CP) BEC condensate-fraction formula ([5]) for d = 3, 2 and 1, provided one assumes s = 1. For 
lower T's, one can argue [35] that a mixed BEC phase containing both 2e- and 2h-CPs becomes more 
stable (i.e., has lower Helmholtz free energy) so that the simple pure-phase formula {8]) is no longer strictly 
valid. Indeed, ([8]) applies at all to the CPs because in the binary boson-fermion gas mixture — for, say, 
a Cooper/BCS model interaction forming the bosonic CPs with a maximum allowed [39] [40] coupling of 
A = 1/2 — only a minuscule fraction (< 0.1%) |41j of the individual fermion charge carriers are paired up 
into CPs, ensuring that a substantial Fermi sea is still present. Such tiny fractions are consistent with some 
very recent far- infrared charge-dynamics measurements [42 43J in LSCO. 



To conclude, we have presented normalized order-parameter data based on penetration-depth and gap mea- 
surements that strongly suggest a linear energy vs. center-of-mass-momentum (dispersion) relation for 
Cooper pairs in various materials that can be viewed as 3D, quasi-2D and quasi- ID superconductors (SCs). 
The linearity is a manifestation of the Fermi sea background in which the pairs propagate, as opposed to the 
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Figure 1: BE condensate- fraction curves for bosons in d = 3, 2, or 1. 



quadratic relation of composite bosons moving in vacuo. It ensures that a BEC picture of SCs is applicable 
over all dimensionalities in which SCs occur. 

In Figure 1 we show BE condensate-fraction curves (in thick) 1 - (T/T c ) d / S for bosons in d = 3, 2, or 
1 assuming dispersion relation (6) for s — 2 and 1, for a pure phase of either 2e- or 2h-CPs, compared to 
empirical data for 3D SCs (Nb/Cu and Sn); for two quasi-2D SCs (Y123 and Bi2212 with T c ~ 93 K and 91 K, 
respectively); and a quasi- ID SC (4- A- wide nanotubes with T c ~ 15 K). The dashed curve labeled d/s = 3/2 
appears in Fig. 6 of Ref. pQ and seems to provide the only adjustable-parameter-free comparison with 
experimental data in that paper. The ordinate axis refers to a universal "normalized order parameter." Data 
for the 3D and 2D SCs refer to penetration depth measurements. Nanotube data are gap A(T) measurements 
giving A(T)/A(0) but are plotted as [A(T)/A(0)] 2 so as to coincide with the 2h-CP condensate fraction 
mo(T)/mo(0) according to the relation A(T) = j \J too(T'), with / a boson- fermion coupling constant (which 
drops out from the normalized order parameter), that follows for 2h-CP condensates from the generalized 
BEC theory of Ref. [36]. The dotted straight line marked d/s = 1/1 strictly corresponds from (7) to T c = 0; 
however, it serves as a lower bound to all curves with d/s = (1 + e)/l > 1 for small but nonzero e, implying 
quasi-lD geometries for which T c > 0. Also shown for reference are the two-fluid model [37] curve 1 — (T/T c ) 4 
and the BCS normalized gap A(T)/A(0) order-parameter curve [35]. 
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